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Abstract. We develop a new variant of the recently introduced stochastic transfer- 
matrix DMRG which we call stochastic light-cone corner-transfer-matrix DMRG 
(LCTMRG). It is a numerical method to compute dynamic properties of one- 
dimensional stochastic processes. As suggested by its name, the LCTMRG is 
a modification of the corner-transfer-matrix DMRG (CTMRG), adjusted by an 
additional causality argument. As an example, two reaction-diffusion models, the 
diffusion-annihilation process and the branch-fusion process, are studied and compared 
to exact data and Monte-Carlo simulations to estimate the capability and accuracy of 
the new method. The number of possible Trotter steps of more than 10^ shows a 
considerable improvement to the old stochastic TMRG algorithm. 
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1. Introduction 

The density-matrix renormalisation group (DMRG), developed by White in 1992 [|l|, is 
one of the most precise numerical algorithms for the investigation of low-dimensional 
strongly correlated systems. Originally the DMRG was introduced to compute the 
ground state and low energy spectrum of a quantum Hamiltonian H. Meanwhile there 
are a number of variants using the basic DMRG idea of numerical renormalisation in 
other physical fields 0. An important progress was made by applying the DMRG 
to the transfer-matrix of 2D classical systems 0, a method known as transfer- matrix 
DMRG (TMRG). This even allows to analyse the thermodynamics of ID quantum 
systems by mapping the partition function to a 2D classical model using a Trotter- 
Suzuki decomposition 0. A highly efficient realisation of a TMRG algorithm, which 
uses corner-transfer-matrices (CTMs), was proposed in |^ and is called corner-transfer- 
matrix DMRG (CTMRG). 
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An upcoming new application field of the DMRG algorithm are ID stochastic 
systems. The dynamics of such models are described by a Master equation which has the 
form of a Schrodinger equation with a "stochastic" Hamiltonian if |]^, ^. The DMRG 
algorithm can be used to compute the stationary limit of the stochastic process, which 
corresponds to the ground state of H. In contrast to quantum systems H is in general 
not hermitian since there is genuinely no detailed balance. Carlon et al first applied 
the stochastic DMRG algorithm to various react ion- diffusion models and discussed in 
detail the influence of non-hermitian operators on the numerics. 

An alternative approach to stochastic models using a TMRG algorithm was 
proposed in |T^. In complete analogy to quantum systems the dynamics of ID stochastic 
models can be mapped to a 2D classical model. Therefore it was quite natural to apply 
the TMRG to the corresponding "stochastic" transfer-matrix. Even though this so- 
called stochastic TMRG is similar to the quantum case in many respects, some important 



differences appear. Enss and Schollwock [|TT| discussed in detail properties of the 
stochastic transfer-matrix focussing on the choice of the density-matrix. Unfortunately, 
the stochastic TMRG shows an unsatisfactory convergence caused by inherent numerical 
problems which are related to the structure of the stochastic transfer-matrix. 

The present work proposes a new approach to analyse the dynamics of stochastic 
problems, a method which we refer to as stochastic light-cone CTMRG (LCTMRG). 
As suggested by its name, the LCTMRG combines ideas of the stochastic TMRG and 
CTMRG algorithms, adjusted by a causality argument which demands a number of 
modifications for an adaption of the CTMRG to stochastic problems. We show that 
the LCTMRG is a considerable improvement of the stochastic TMRG with respect to 
numerical stability and performance. 

2. Stochastic Models 

Stochastic models have gained a large interest in statistical physics. They are used 
not only in physical but many interdisciplinary research fields to describe processes far 
away from thermal equilibrium [M . The bandwidth of applications reaches from the 



description of social behaviour and biological processes to traffic flow (see e.g. |13, p!4|). 
Typically, stochastic models start from an initial state which evolves in time according 
to (local) probabilistic rules. 

In the present work we focus on one-dimensional stochastic problems. We consider 
a chain of length L, where each site Si can either be occupied by a particle (sj = 1) or 
empty (sj = 0). In stochastic physics one is interested in the dynamic evolution of a 
probability distribution P{t) of states. P{t) can be denoted as a vector 

\P{t)) = Y.Ps{t)\s) (1) 

s&S 

where Ps{t) is the probability of finding the chain in the configuration 

s = (si,S2,...,Si)G5 = {0,ir^. (2) 
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Depending on the type of dynamics, \P{t)) can evolve in continuous or discrete 
time. Assuming continuous dynamics, stochastic processes can be described by a Master 
equation 

d,\P{t)) = -H\P{t)) . (3) 

H is called "stochastic Hamiltonian" [|^, § because has the form of a Schrodinger 
equation in imaginary time. The matrix elements are given by 

{s\H\s) = -w{s s) + 6s,s ^ w{s s') (4) 

s'es 

where w{s —>■ s) denote the probabilistic rates of the transition s ^ s. Although the 
Master equation suggests a close analogy to quantum systems, as an important 
difference the stochastic Hamiltonian H is in general not hermitian. A formal solution 
of is given by 

|P(t))=e-*-^|P(0)) (5) 

where \P{0)) denotes the initial probability distribution at t = 0. Obviously, the 
stationary limit \P{oo)) is a (right) eigenvector of H with eigenvalue 0. 

Stochastic models can show a rich phase diagram and interesting critical phenomena 
The simplest situation is an absorbing phase transition [|I2| into an empty state 



|00- ■ - 0). Introducing the occupation number operator nj of site j, the average local 
density of particles 

n{t) = {l\nj\P{t)) with (1| := ^ {s\ (6) 

ses 

(with arbitrary j due to translational invariance) is an order parameter that distinguishes 
the phases: all particles can either vanish (n(oo) = 0) and the system falls into the 
(absorbing) state |P(oo)) = |00---0), or a certain number of particles stay "active" 
(n(oo) 7^ 0). In the critical region n{t) evolves according to a power law whereas a 
non-critical behaviour is characterized by an exponential decay: 

{e"*"^"^ non-critical 
(7) 
t " critical . 

Like quantum systems, stochastic processes show universal behaviour at criticality. 
The most prominent universality class, which is typical for phase transitions to an 
absorbing state, is the directed percolation class (DP) |T^. But in general non- 
equilibrium phase transitions are by far not so well understood as those in equilibrium 
physics. 

3. Stochastic Transfer-Matrix DMRG 

3.1. The Stochastic Transfer- Matrix 

We focus on the calculation of the dynamic evolution of the local density 

n(t) = (l|n-e-*-^|P(0)) (8) 
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of a stochastic process in the thermodynamic hmit L — » oo and assume that H consists 
of local "stochastic interactions" 



(9) 



As in the conventional quantum TMRG algorithm [Q, the stochastic system is first 
mapped to a 2D statistical model by using a Trotter-Suzuki decomposition of (||). 
The resulting classical 2D lattice |T^, |Tl| is shown in figure |I] (a) which local plaquette 
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Figure 1. (a) Trotter-Suzuki decomposition of 2 At time steps. The resulting 2D 
lattice consists of local plaquette interactions r and is infinitely extended in space 
direction. The dimension of the time direction is finite and the boundary conditions 
are fixed by (l| and |P(0)). (b) Reduction of the 2D lattice to a triangle structure. 
All other plaquettes trivialize due to (p^). 



interactions are given by 

h rj 

Mnr. = (^2r2|e-^*'^|/iri) = ^ with /,,r, e{0,l}. (10) 

h ri 

Thus, the spatial dimension L of the stochastic process is expanded by a virtual Trotter 
dimension M = t/At which corresponds to the time direction and is split into (discrete) 
steps of size At. At has to be chosen sufficiently small to obtain a good approximation 
of n{t). Formally, the Trotter decomposition becomes exact for At 0. As we measure 
the local density n{t) at finite time t, but in the thermodynamic limit L ^ oo of 
the stochastic chain, the space dimension of the 2D lattice is infinite, whereas the 
Trotter dimension is finite. Note that in contrast to the quantum TMRG the boundary 
conditions are fixed in Trotter direction and given by the vectors (1| and |P(0)), cf. 
figure (a). 

In complete analogy to the quantum case one can apply a TMRG algorithm to 
the 2D lattice fT^, 0. Using column transfer-matrices, shown here pictorially for the 
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T 



and T{n) 



J 



the local density n(t) for the thermodynamic limit L oo can be calculated by 

n{t) = {il;L\T{n)\i;R) . (12) 

IV'-R/l) labels the leading right/left eigenvector of T, having the eigenvalue 1 |Tl|]. 
The stochastic TMRG algorithm is used to compute iV'iJ/i) and T{n) for successively 
increasing Trotter numbers M. Unfortunately, various computations show numerical 
problems that limit M ~ 10^ [p!0| , |11|, which is far from enough to compete with other 
methods like Monte-Carlo simulations (MCS). 

Here we propose a different TMRG algorithm based on corner-transfer-matrices 
(CTM). Such a corner-transfer-matrix DMRG algorithm (CTMRG) is known to be 
numerical more stable and faster than TMRG 0. Before we explicitly construct these 
CTMs for the stochastic case, we discuss some physical properties of the 2D lattice 
relevant for the development of the new algorithm. 

Due to probability conservation we have (l|e~^* '^|s) = 1 for any state \s). Thus, r 
trivializes by summing out the "future" indices, i.e. 

T.I2 Sra 



V/i,ri: 



r 



(13) 



hr2 



The effect of ( P^ to the 2D lattice is discussed in detail in |T^ . For the computation of 
n{t) it is found that a huge number of plaquette interactions can be omitted, because 
they "trivialize". The remaining non-trivial plaquettes form a 2D lattice of finite 
dimension which is shown in figure |I| (b). The trivialisation process can easily be 
understood by a causality argument: only a "light-cone" of plaquette interactions can 
influence the site where the local density is measured. 

We now construct a CTMRG algorithm which genuinely fits to the triangle structure 
of the 2D lattice. As shown in figure H, four cuts are set to separate the lattice into four 
parts. The cuts are somewhat native to our model, because they form the boundaries 
of the "future" and "past light-cone" of the center point of the triangle. The four parts 
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Figure 2. Construction of corner-transfer-matrices, (a) The lattice is split into four 
parts, (b) Schematic plot of the corner-transfer-matrices C„, (7^,(7;, and Cr evolving 
from (a). 



are interpreted as CTMs whereby and label block-spins. We next show how 
these CTMs can be treated within a CTMRG algorithm analogous to e.g. 0. However, 
a number of modifications are necessary to adapt the CTMRG to the light-cone of 
plaquettes. Hence, we call this CTMRG variant hght-cone CTMRG (LCTMRG). 



3.2. The Light-Cone-CTMRG Algorithm 

In a CTMRG algorithm the CTMs are enlarged sequentially by adding transfer-matrices 
(TMs) to each cut. We first define the TMs in a pictorial way 



Id 



rd 



lu 





(16) 



(17) 



The bullets represent single spin sites s, s and the spins rig, fig marked with bars will 
become renormalized block-spins in the DMRG algorithm. Exemplarily, ([T6| ) and (p!7|) 
show TMs which are used to enlarge the triangle lattice of figure ^ (a). 

Figure |^ demonstrates graphically how the TMs are used to enlarge the CTMs, 
whereby CTMs and TMs are "jointed" by summing out the adjacent indices (similar 
to a matrix multiplication). Due to the exotic geometry of the lattice we have to 
distinguish between an upper and lower extension step depending on whether Cu or Cd 
should be enlarged. In our LCTMRG algorithm both extension steps are implemented 
alternately. That way all CTMs grow step by step with the crossing point of the cuts 
always situated in the center of the triangle. After each extension step the CTMs have 
to be renormalized by a density-matrix projection, cf. section |3.3|. 
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Figure 3. Extension of tlie CTMs by adding diagonal TMs. In (a) the lower part of 
the triangle is enlarged whereas in (b) the extension is performed to the upper part. In 
the LCTMRG algorithm the upper and lower extension are done alternately, so that 
all CTMs grow step by step. 



The local density n{t) can be obtained using Cd, and T^d 



E s 




It is important to notice that n{t) is computed in the center of the triangle-lattice. 
Here, influences of boundary effects are expected to be smallest. In terms of a DMRG 
algorithm, the CTMs C;, Cr and Cu act as an "environment" of the "system" Cd- 



3.3. The Choice of the Density-Matrix 

The key problem is to find a reasonable density-matrix projection for the renormalization 
of the CTMs after each extension step. We exemplify the construction of the density- 
matrix by looking at figure |^ (a). Here, one block-spin and two spins of Cd, Ci and Cr 
have to be renormalized into one block-spin. The construction of the optimal density- 
matrix projection is now discussed in detail. 
First, we define four vectors: 
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RJus 



■S2 



'Si 



3l 



51, S2 



(19) 
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L)ns,si,S2 




(20) 



The block-spin fis belongs to the environment, si, S2 and to the system block. Note 
that these vectors approximate the left and right eigenvector of the leading eigenvalue 
of diagonal TMs 



■R 




and Tl= ^ (21) 
which have a different shape compared to the stochastic TMs used in the stochastic 




TMRG [jTO|, |TT|]. As and V'l = -R) have not to be computed by any expensive 
diagonalization routine like in TMRG, the LCTMRG algorithm is much faster. 

■0/? and ipf^ are used to construct a reduced density-matrix for each of the two cuts. 
The most generic choice would be a symmetric density-matrix 

f^^ = trnAmm + \rR)m) m) 



which was also used in [^, |T0|, [TT|. Here tr^j^ denotes the partial trace over n^. pL^^ 
produces a reduced system block basis which optimally approximates and 0. 
However, one can easily proof that ip'f^^ is trivially given by 

= 1 for alln„Si,S2,n, (23) 

which follows directly from the trivialization process (f^). Obviously is not very 



useful for constructing a density-matrix, because 

trsj^i?)(^i?| = with (l.)„„.,,,, = 1 (24) 

reduces to a trivial projector which does not correlate system and environment block. 
Hence, we omitted and tested the density-matrix 

Pi''=tr.j^2)(^2| (25) 

which led to much better results (cf. section |4^ ). An asymmetric choice 

pf=^^nArR){rL\ (26) 
of the density- matrix performs worst. As 

is independent of ?t,^,s'i,S2, the density-matrix px has rank one and represents a pure 
projector. 

[21 

A physical explanation for the choice of px can be given in terms of the light-cone 
picture of section ^7\\ . The trivial vector \ip^ is a superposition of all feasible states 
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which means that at each cut no further information about the "future" is available. 
Not surprisingly, the system and environmental part of the cut are uncorrelated which 
is expressed by a trivial density-matrix projection of tifis \'^R){''Pfi\- Only the physics of 
the past, the information of which is carried by IV'l)' correlate system and environment 
indices. 



3.4- Some Technical Aspects 



We briefly discuss some implementation details of the LCTMRG algorithm. 
The first TMRG step should start with the following configuration 
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of transfer-matrices. This corresponds to the time t = 1.5 ■ At and is the simplest 
construction of initial CTMs and TMs. The first extension steps (cf. figure ^ are 
performed without renormalisation until the dimension of the CTMs exceeds the number 
m of DMRG states. 

As the transfer- matrix is only used for the construction of ipf^, which is not 
needed for computing px , Cu can be omitted completely. If additionally the local 
Hamiltonian is parity invariant, i.e. = the local transfer-matrix r 

becomes symmetric. Hence, only the CTMs Ci,Cd and TMs Tid,Tiu have to be stored. 
Cr,Trd and Tru can be reconstructed by mirroring Ci^Ti^ and T;„. 

In order to avoid floating point overflows of the algorithm, it is recommended to 
rescale all CTMs and TMs 



Cx 



Cx 



T 



T 



(28) 



liar \\TA 

where || ■ || is some norm. Note that these rescaling factors have to be considered in the 
computation of n{t). 

All computations were done on Sun Workstations (Ultra Sparc III, 900 MHz). 



Compared to the old stochastic TMRG [0, |TT| the LCTMRG algorithm is tremendously 
more efficient. Furthermore, as most parts of the program consist of matrix multiplica- 
tions of CTMs and TMs, the LCTMRG algorithm can easily be parallized. The CPU 
time needed for one Trotter step ranges from a few milliseconds for m = 32 states up 
to a couple of seconds for m = 400 states. The consumption of computer memory is 
modest as well, e.g. 10 MB for m = 32 up to 200 MB for m = 400. 
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4. Application 

4.1. The Model 

The numerical studies of the present work focus on react ion- diffusion processes (RDP) 
which are used to model various chemical reactions. We consider a simple RDP with 
one type of particle A which exhibits the following reactions 



(29) 



This RDP can be expressed by a stochastic Hamiltonian H = ^ . /ij^j+i with local 
interactions 

f -S ~6 -2a \ 

D + S + P -D -7 
-D D + 5 + P -7 
\0 -f3 -(3 2(a + 7) / 

li^i+i is parity invariant and the local transfer-matrix r becomes symmetric. Therefore 



diffusion: 


AH) 
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(with rate D) 


annihilation: 


AA 


^ 00 


(with rate 2a) 


coagulation: 


AA 


^ U,A(D 


(with rate 7) 


death: 


U^AdS 


^ 00 


(with rate 6) 


decoagulation: 


0A^0 


^ AA 


(with rate P) 



h 



(30) 



the LCTMRG algorithm simplifies, cf. section p^ . 

As an example, we apply the new LCTMRG algorithm to two RDPs, the diffusion- 
annihilation process and the branch-fusion process. These models have also been chosen 
by Carlon et al |^ to demonstrate the efficiency of the stochastic DMRG algorithm. 

The diffusion-annihilation process (DAP) 

2D = a, /5 = 7 = 5 = (31) 

is exactly solvable Jl^. For an unbiased initial probability distribution the dynamic 
evolution of the local density is given by 

nit) = \ (loiADt) + h{ADt)) e"^^* (32) 

where /q, h are modified Bessel functions. Thus we can use analytical results to check 
the numerical precision of the TMRG data. Note that the DAP is critical for all D with 
an asymptotic behaviour n{t) ~ t^^/^. 
The branch-fusion process (BFP) 

D = 2a = -f = 5 =: 1-p, (3 =: p (33) 

is a simple one parameter model which exhibits a non-equilibrium phase transition 
from an active to an absorbing phase. The BFP is not exactly solvable, but Monte- 
Carlo-Simulations and stochastic DMRG computations ||^ are available. The critical 
behaviour of the BFP falls into the DP universality class where precise data for the 
critical exponents have been calculated by series expansions [TH . 



Using the stochastic LCTMRG for both models can demonstrate, whether the 
method is capable to produce reliable results for 
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• critical and non-critical systems 

• systems at phase transition points. 

The next subsections present numerical results for the two processes. If not stated 
differently, all computations were performed with At = 0.05. 



4-2. The Diffusion- Annihilation Process 




Figure 4. The left figure compares exact data for the DAP with LCTMRG 
computations by showing the absolute error e = |nLCTMR,G(i) ~"'exact(i)| for D = 0.05, 
keeping m = 32, 64, f 00 states. The inset plots LCTMRG data for D = 0.05 and 
TO — 100 in a double-logarithmic plot and shows the algebraic decay. The right figure 
depicts LCTMRG data for D = 0.5, m = 32 by using different density matrices p'^' 
and 



Figure ^ (left) compares LCTMRG calculations with exact data for D = 0.05, 
keeping various numbers of states m. Up to more than M ~ 10^ Trotter steps we 
obtain highly precise data with a deviation of less than 10~^ from the exact results. 
The inset of figure § (left) plots the LCTMRG data in a double-logarithmic plot, which 
shows that n{t) falls of algebraically. 

The high number of Trotter steps M is a considerable improvement to the old 
stochastic TMRG algorithm [|l^, ^ by at least three orders. Even though the DAP is 
critical, we observe an extremely stable convergence of the LCTMRG algorithm. 

Figure § (right) plots numerical data for different density-matrices p^^^ and p'^', 
cf. section |3]^. In all our calculations we observe highly instable numerics, if the 
conventional density-matrix p'^1 is used. In the example of figure ^ the convergence 
of the algorithm breaks down after M ~ 10^ Trotter steps, while M ~ 10^ is possible 
for p[^]. Thus the arguments given in section |3.3| can be confirmed numerically: p[y is 
not an adequate density-matrix for the stochastic LCTMRG. 
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4-.3. The Branch-Fusion Process 

In this section we focus on the critical phase transition of the BFP at pc = 0.84036(1) 
Figure H (left) compares numerical data computed by the LCTMRG algorithm with 
conventional Monte-Carlo simulations. For p sufficiently far away from criticality, we 
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500 1000 1500 

time M=t/At 



Figure 5. The left figure siiows the dynamic evolution of the order parameter n{t) of 
the BFP for p = 0.9 . . . 0.8 (in steps of 0.01) keeping m = 200 states. For comparison 
Monte-Carlo simulations (o) are plotted. The right figure plots the logarithmic 
derivative C{t) which converges to the critical exponent a. The literature value is 
plotted by a dashed line. The inset shows C{l/t) which is used to interpolate the 
critical exponent a. 



observe a convergence up to more than 10^ Trotter steps. In the vicinity of the critical 
point p ~ Pc the convergence becomes worse. Figure ^ (right) plots the logarithmic 
derivative 

^ logn(t + At) -logn(t) 
^ ^ logAt ^ ' 

at the critical point p ^ Pc for m = 400 states, which is very sensitive to numerical 
errors. Up to more than 10'^ Trotter steps the numerics are extremely precise, and one 
can verify that n{t) switches to an algebraic behaviour. It is also possible to determine 
the critical exponent a by extrapolating C(t ^ oo), cf. inset of figure | (right). Thereby, 
we were able to compute a up to a precision of less than 0.1%: 

a ^ 0.1600(5) (literature a^t = 0.159464(5), [18]). (35) 

The question arises why the convergence of the LCTMRG at p ~ is two orders 
less than in the DAP process, yet both models behave critical. We checked various 
numerical aspects to determine the origin of the worse convergence. 

In the BFP it is conspicuous that the quality of the results at p ~ Pc strongly 
depends of the number of states m that are retained within the LCTMRG algorithm. 
This is demonstrated in figure ^ (left) which plots numerical computations for various 
m at p ~ Pc- However, it is surprising that so many states m are needed although there 
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Figure 6. The left figure plots the dynamic evolution of the order parameter n(t) of 
the BFP at criticality p ^ Pc for various number of kept states m = 32 . . . 400. A MCS 
is plotted for comparison by a dashed line. On the right side one can see the spectrum 
Pi of the density-matrix p for Trotter steps M = 10, 50, 100, 1000 at criticality p ~ Pc- 




Figure 7. Both plots show the order parameter n{t) of the BFP at criticality p ^ Pc, 
keeping m — 6A states. On the left side At = 0.001 . . .0.5 is varied. The right figure 
compares data calculated with various floating point precisions (i.e. 64, 96, 128 and 
160 bits). The plot shows the differences of numerical data for each precision compared 
to the next higher one. 



is a very strong fall off of the density-matrix eigenvalues pi, cf. figure ^ (right). 

As another possibility we check the influence of the size of the Trotter steps At on 
the numerics. The curves of figure ^ (left) belong to various At, but are rescaled to 
At = 0.05 for comparison. Even though finer Trotter decompositions increase the total 
number of convergent Trotter steps, one can not improve the accuracy of the data with 
respect to the absolute time t. If on the other hand At becomes too large, the Trotter 
decomposition itself gets worse and is then responsible for unsatisfactory numerical data. 
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To estimate the effect of numerical errors caused by floating point inaccuracies we 
implemented the diagonalization routine for the density-matrix alternatively by using 
higher mantissa bits. This was technically realized by using the GMP library [1^] which 



allows an arbitrary number of mantissa bits. As shown in figure ^ (right), only a marginal 
effect on the numerics is observed. 

Overall, it remains an open question what exactly is the limiting factor for the 
worse convergence at the phase transition point p ~ in the BFP. To exclude model 
specific reasons, we also checked other RDPs, e.g. the contact process. Qualitatively, 
the same limited convergence near the critical phase transition point is observed. 

5. Conclusions and Outlook 

In the present work we proposed a new variant of the stochastic TMRG by using 
corner-transfer-matrices which we call stochastic light-cone CTMRG (LCTMRG). The 
LCTMRG algorithm fits genuinely to the specific structure of the triangle classical 
lattice which evolves from the Trotter- Suzuki decomposition of the stochastic model. 

We tested the new algorithm by comparing LCTMRG data to exact results and 
Monte-Carlo simulations of two different reaction-diffusion models. We obtained highly 
precise numerical results (e ~ 10~^) up to M ~ 10^ Trotter steps, even if the model 
behaves critical as the diffusion-annihilation process. Compared to the old stochastic 
TMRG algorithm ITO, |1T| with M ~ 10^, this is an enormous increase of the number 



of reachable Trotter steps of three orders. An important observation is that inherent 
numerical problems of the old stochastic TMRG algorithm do obviously not appear in 
our new approach. In the vicinity of a critical phase transition point, exemplified by 
the branch-fusion model, the convergence gets worse, but is nevertheless sufficient to 
determine precise results for critical exponents. 

Since within critical phases a much better convergence has been observed, it 
remains an open question what exactly causes the reduced convergence at critical 
phase boundaries, which presumably does not originate from purely numerical reasons. 
Therefore our future research is concentrated on further modifications and improvements 
of the LCTMRG algorithm, e.g. the implementation of a finite size algorithm. 

Overall the numerical investigations show that the new LCTMRG algorithm is a 
considerable step towards a general and very efficient method for ID stochastic problems. 
Compared to the traditional approach using Monte-Carlo simulations, there are two 
fundamental advantages of the LCTMRG: 

• The LCTMRG is not a simulation technique. There is no need of taking random 
numbers and sample averages. The LCTMRG is a numerical renormalisation group 
based on the quantum formalism for stochastic models where averages are directly 
accessible. 

• The algorithm describes the exact thermodynamic limit L — oo of the stochastic 
model. Note that here we even have to deal with a finite classical 2D system 
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only, due to the simplification from the "light-cone decoupling" . Thus, there are in 
principle no finite-size effects like in MCS or stochastic DMRG. 

Even if the number of possible time steps, in particular at phase transition points, can 
not compete with MCS up to now, we believe that stochastic TMRG algorithms can 
be an extremely valuable tool for studying ID stochastic systems. Finally we mention, 
that - as in the case of the CTMRG - a generalisation of the LCTMRG to more 
than one dimension is also imaginable. 
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